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1. Introduction 

Protoneutron star is a hot, lepton-rich neutron object which is formed in Type-II 
supernovae explosion as the final stage of a star more massive than about 8 M & - 
25 M Q |1J. The evolution of a protoneutron star depends on two factors: the relativistic 
equation of state of the stellar matter and neutrino interactions which provide to the 
phenomena of neutrino trapping in a neutron star matter . Neutrino interactions play 
a crucial role during the creation of the protoneutron star. At the beginning, this star 
is very hot and lepton-rich, after a typical time of several tens of seconds, the star 
becomes cold, deleptonized object and the neutron star is formed. During the collapse 
enormous energy of the order of ~ 10 53 ergs is released just through neutrinos. The 
released energy is equal to the gravitational binding energy of a newly formed neutron 
star. In this paper the mass-radius relation of the protoneutron star is constructed for 
different temperature cases. A particularly interesting situation results in the occurrence 
of unstable areas which allows distinctively different evolutionary tracks. In the new- 
born protoneutron star (less than several seconds) neutrinos are trapped locally in the 
dense stellar matter at densities greater than 10 13 g/cm 3 , forming an ultrarelativistic 
and degenerate Fermi gas. This trapped neutrinos have an influence on the equation of 

state HBBgg. 

The outline of this paper is as follows. 

In Sect.l the general properties of nonrotating and rotating protoneutron stars are 
calculated. In Sect. 2 the employed equation of state is obtained in the approach of 
relativistic mean-field approximation. This approach imply the nucleons interactions 
through the exchange of meson fields. Thus the model considered here comprises: 
electrons, muons, neutrinos and scalar, vector-isoscalar and vector-isovector mesons. 
Consequently the pressure and energy density are characterized by contributions coming 
from these components. In Sect. 3 the influence of rotation on the protoneutron stars 
parameters such as: masses, radii and moment of inertia are calculated. The obtained 
results for the given form of the equation of state followed by a discussion of the 
implications of these results are presented in Sect. 4 and Sect. 5. 

In this article the signature { — , +, +, +} is used, in which the Dirac's matrix is defined 
by 

and 

and h — c — 1. In case of the fermion fields it is more convenient to use the reper field 
defined as follows g^ v = e^e^n^ where i] ab is the flat Minkowski space-time matrix. 
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2. The theoretical model 

This paper presents a basic model of a protoneutron star in the RMF approximation 
H 0; PI P; FH- The relativistic approach to high-density nuclear matter in this approach 
was proposed by Walecka. Since then, this theory has been applied successfully to many 
subjects in nuclear many body problems. Especially the Walecka model (QHD) and its 
nonlinear extensions have been quite successfully and widely used for the description 
of hadronic matter and finite nuclei. Increasing interest in the neutron star matter at 
finite temperature has been observed recently in relation to the problems of hot neutron 
stars and protoneutron stars and their evolution in particular. Theories concerning 
protoneutron stars are being discussed in works by Prakash et al. |], ||, |], [II], [TJ|, |I3" 



Glendenning |14| has studied the properties of neutron star in the framework of nuclear 



relativistic field theory. In the field theoretical approach the nuclear matter is described 
as the baryons interaction through mesons a, uj and p exchange. We define an action of 
the protoneutron star through the Lagrangian density £ 

S = J al 4 Xy/^C. (1) 

The Lagrangian density function of this theory can be represented as the sum 

£ = C B + C L + C M + C G , (2) 

where Cl, C-m-, £-g describe the baryonic, leptonic, mesonic and gravitional sector, 
respectively. The fermion fields are composed of neutrons, protons, muons, electrons 
and neutrinos 



The baryonic sector of the Lagrangian density function is given by 

C B = ilfl-fD^ - Tp(M - g a cp)iP, (4) 
where _D M is the covariant derivative 

D„ = d^ + l;ig P p*o- a + ig^u^. (5) 
The leptonic part of the Lagrangian (Q) is define as 

C L =iJ2L fl »d„L f -J29f(L f He Rf + h.c), (6) 
/ / 
where H is Higgs field and this field has only the residual form 

°) 

the value of V = 250 GeV comes from the electroweak interaction scale. In the mesonic 
sector the Lagrangian density function is given by 
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9a 


9w 


9 P 


92 


MeV 


93 


10.0289 


12.6139 


4.6322 


1427.18 


0.6183 


M a [MeV] 


M^\MeV\ 


M p [MeV] 


M 


MeV] 


c 3 


511.198 


783 


770 


938 


71.3075 



Table 1. The parameters set of the model [fig 



where <p is the scalar field, F pv is the electromagnetic stress tensor 



flV 



d p A u — d v Ap 



and Wfiv and R a are vector mesons field strength given by 



W, 



fJ,V 



pa 

K fJ,U 



%Pl ~ d ^P p + 9 P £ahcP»Pl 

M w and M p are the meson u and p masses. The potential function U((p) possesses a 
polynomial form 

U(cp) = l -M 2 y + lg 2 <p 3 + Igw* 



introduced by Boguta and Bodmer [15] in order to get a correct value of the 



compressibility K of nuclear matter at saturation density. The model described by 
the Lagrangian function (g) with such an ansatz for the potential U(ip) guarantee a 
very good description of bulk nuclear matter properties for different parameter sets. It 
is not possible for neutron star matter to be purely neutron one. As it was stated above 
the fermion fields are composed of protons, neutrons, electrons, muons and neutrinos 
because we are dealing here with the electrically neutral neutron star matter being in 
/^-equilibrium. Such a matter possesses a highly asymmetric character caused by the 
presence of small amounts of protons and electrons. Finally, the standard gravitational 
sector has the form 

1 



C 



G 



2k 



-R 



(7) 



where k = 8ttG/c a . The parameters entering the Lagrangian function are the coupling 
constants g w , g p and g a and self- interacting coupling constants g 2 , g%- The parameters 
employed in this model are collected in the Table |l|. In our calculations, we used the 
TM1 |16 parameter set, which has a capability to reproduce the known results of finite 
nuclei as well as normal nuclear matter. The equations of motion obtained from the 
Lagrangian function for individual fields are the Dirac equations for nucleon fields tp 



ij^D^i/j — mpip = 
with mp being the effective nucleon mass 

m F = MS = M - g a ip 
and Klein-Gordon equations with source terms for meson fields (p, uj^, p p 

U(p = Ml if + g 2 y 2 + g 3 p 3 - g a ipi; 



(8) 



(9) 



(10) 
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= M%u" + csK^K + 9u,r B (11) 
D^R^ a = M 2 pP va + g p J au . (12) 
Sources that appear in the equations of motion are the baryon current 

r B = ^ (is) 

and existing only in the asymmetric matter the isospin current 

J av = W\<j a ^- (14) 

The baryon and isospin charges are defined by the zero component of the adequate 
currents, so 

Qb = J d 3 xn B = J d 3 xJ° = J d 3 xip + ip = J d 3 x(ip+ if; p + ^+ ip n ) (15) 

and 

Q 3 = J d 3 x^°^a 3 iP = 1 J d 3 x(^ p - ^n). (16) 

The form of the energy-momentum tensor that results from the variational principle 
is given by: 

T n dC M dC F 
^ ~ dg» u M de au ^ matt ' 1 ' 

The aim of this paper is to achieve the equation of state of the protoneutron star 
matter at finite temperature. Our calculations are based on the variational method 
incorporating the Feynman-Bogoliubov inequality. They were presented in details in 
paper [10]. The total pressure of the protoneutron star appearing in T pv is the sum of 
fermion and meson parts P = Pf + Pm-, 

p = p ^KM^ + 1^ + ^{S) 2q1 - u(v) - (18) 

Whereas the total energy density coming from fermion and meson contributions e = 
Bp + Em can be written as 

£ = £ F-l^w) ^-\c^ + 9^Q B + \gl[^ QI + U&). (19) 

The fermion pressure and the energy density in the case of massive particles are defined 

as 



1 r°° k^ / 

V ' / dkk 2 ; (fj+J. 

i 3tt 2 Jo Jk* + mh K 



] poo , 

^ = E-/ dkk 2 J¥+ri F ~(f i + f i ) 

; TT Z JO v 



where fa and f\ are the Dirac fermion distribution functions. The pressure (-Pf) and the 
energy density (ep) can be written with the use of the functions 0, and Xi as 

P F ,i = PoM x i)> (2°) 
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£F,i — £ 0,iXi{ x i) 

where 

n m F)i c 2 

A is the Compton wavelength. The dimensionless functions fa and Xi are stated in the 
following way 

fa(z,t) = — / A fi(x,z,t) + f i (x,z,t)}= (21) 

3tH Jo x 2 + 



37T 2 Ay' 



Xi ( z ,t) = —j o x 2 dxjx^6 2 {f i (x,z,t) + f l (x,z,t)} = (22) 

= OHbO", + m 2 H 3 (r, yi)) 
with the distribution functions 

*) = 7 ; r \ , (23) 

7*(*,«,*) = - r - F =^ • (24) 

e (v^+«)/* + 1 

In those relations = M5j, x = kp/M, z = /x/M, r = z/t = fx/hsT, yi = 

Si/t = mF,i/kBT and t = ksT/M are the effective mass, the dimensionless Fermion 
momentum, chemical potential and temperature, respectively. A^ = 1/ksT is the 
thermal wavelength (5j = 5 for i = p,n end <5 e = m e /M where M is the neutron 
mass). The thermodynamic properties of fermions can be evaluated using the well 
know functions H n and G n given by 

jf.= i r-p*_{ /+ 7} (25) 



37r2 Jo ^Jk 2 + m 2 F 
G n = - 2 j o k n ~ l dk{f-f) (26) 



The two terms in (|25|j26|) correspond to the contributions of particles and antiparticles, 
respectively. Our first step is to calculate the pressure and energy density for nucleons 
which correspond to the nonrelativistic case. The relativistic case, corresponding to the 
high temperature approximation, is utilized in order to calculate the pressure coming 
from electrons and massive neutrinos. In this approach the function H 5 and H 3 are 
given by the following relations: 

H 3 (r, y) = l -y\-l + 2r 2 + 2T) + ^y 2 ln(yM (27) 



and 
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Figure 1. The equation of state for different temperature cases. 



In the case of massless neutrino (m u = 0) the pressure P and the energy density e can 
be expressed analogously to the massive fermion relation (^). In this situation 

Pou = kgT 4 

and function is given by 

1 1 f°° { 1 11 

Mz,t) = -Xo(z,t) = — I x 3 dx\ ,,. t ■ - + - ! -— Tt — . (29) 



3 auv , j 37[2 j o { e (x-z)/t + 1 e (x+z)/t + 1 

The function cf) can be written with the use of the function H n defined as 
H n (r, y = 0) = f X ^d X {^^j + 3:^} 

In agreement with previous substitution the variable r = z/t = fi/kgT and z = fi/M. 
In this moment the case y = is considered. Having integrated these functions (^) one 
can obtain the results in which the solution is expressed by the polylogaritmic function 

Li n (z) 

fa{z, t) = -— 2 {L U {-z~ zlt ) + Li A {-e z / 1 )). (30) 

In the case when the temperature (T = 0) the integral fa and Xi gi ye as a solution the 
well-known Shapiro resalt [ I7| : 

4>{x) = ^ ^xVT+x* ^ - lj + ln(x + vT+^)| , 

x (x) = — {xVl + x 2 (l + 2a; 2 ) - ln(x + Vl + £ 2 )} . 
The total function is sum of the partial components 

= n + (f) p + (j) e + 
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The protoneutron star can be characterized by two parameters, the lepton fraction Yl 
and the entropy per nucleon S. The fermion number density and the entropy per nucleon 
are giving by the relations 



71; 



\ roo , . 

— J x 2 dx{fi(x,z,t) - fi{x,z,t)^ , (31) 



n B 

while the fermions entropy density s is defined as 

\ roo 

s = V— / dxx 2 [-fi(x,z,t)lnfi(x,z,t) - (1 - fi(x,z,t))ln(l - fi(x,z,t)) 

i 7T Z JO 

-7i(x,z,t) hxj^x.z.t) - (1 -7i(x,^,t))ln(l -Ji(x,z,t))] 

|18| with the use of the previously defined Dirac distribution functions fa and f { . 
Considering the case of massless neutrinos the equation ( pi]) can be solved analytically 
and the final form of the neutrino number density is giving by 



n 



xH 
2t? 



ln((l + e ( ^ )/t )/(l + e^/*)) + 



— 2 £ Li 2 (-e^l t ) + - 2 J2 Li s (-e^>/*) 

7T — 7T Z — 

?=(-,+) ?=(-,+) 
The obtained result allows us to determine the lepton fraction Yl, defined as 

Y L = , (32) 

n B 

where n B , n e and n u are the nucleon, electron and neutrino number densities, 
respectively. The proton fraction can be defined as 

Y P = ^, (33) 
n B 

where n p is proton number density The difference between an ordinary neutron star 
and a protoneutron star is caused mainly by the high lepton number of the protoneutron 
star matter. This is because a star core is opaque to neutrinos. The lepton number is 
approximately constant after core bounce Yl — 0.4 for the densities above 10 12 g/cm 3 
|Ij] - 10 13 g/cm 3 H where the neutrinos are trapped. At the time the entropy per baryon 
is approximately constant S ~ 1. Below this density the neutrinos are free and they 
can escape. After the leak out of neutrinos their chemical potential vanishes. The 
situation, when the neutrinos are trapped, is characteristic to very initial stage of a 
protoneutron star existence. The matter is assumed to be composed of nucleons and 
leptons. Thus we are dealing here with electrically neutral matter being in j3 equilibrium 
which can be expressed as a relation between the chemical potentials of the protoneutron 
star components. 

fin + He = /V + /V (34) 
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Figure 2. The lepton number Yl dependence on the central density p c for the different 
temperatures and for the constant proton fraction Yp = 0.1776. 




Figure 3. The entropy S dependence on the central density p c for the different 
temperatures and Yp = 0.1776. 



This equation can be written in terms of neutron x n and proton x p Fermi momentum 



i 



Yp \T 



Si + xl + J SI + xl = M +xl + \I51 + x\ 



for electrons and neutrinos (i = e, v) b% = m e /M,0 whereas for nucleons 5 = 5 n = S p . 
Fermi momentum can be effectively only a function of the neutron Fermi momentum x n 
and the Yp parameter, which measures the proton admixture in the neutron star. The 
protoneutron star matter in general is more symmetric than neutron star matter. 
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Figure 4. The pressure P versus the radius r for different components of pressure. 
The temperature is fixed T = 50 MeV, p c = 1.8 10 15 g/cm 3 . 
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Figure 5. The lepton number Yj, dependence on the radius r for the temperature 
T = 50 MeV and p c = 1.8 10 15 .g/cm 3 . 



3. Rotating protoneutron star 

The spacetime outside a rotating protoneutron star is much more complicated than 
the metric outside a non-rotating star thus it seems to be interesting to investigate 
the properties of rotating protoneutron stars. In general the metric of a stationary, 
axisymmetric, asymptotically flat spacetime has the form 

ds 2 = g^dx»dx u = -e 2u c 2 dt 2 + e 2 ^(d<p - ucdt) 2 + e 2 ^d6 2 + e 2X dr 2 , (35) 

with g^y being the metric tensor. The metric potential functions v, ?/>, fi and A and the 
angular velocity to of the stellar fluid in the local inertial frame are functions of the radial 



coordinate r and the polar angel 9 |22, [23[ . In order to compute the structure of rapidly 
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Figure 6. The entropy S dependence on the radius r for the temperature T = 50 
MeV and p c = 1.8 10 15 g/cm 3 . 
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Figure 7. The mass-radius diagram for the protoneutron star for constant Y& = 0.4. 

rotating fluid body the numerical method developed by Butterworth and Ipser was 
introduced. Besides this exact numerical treatment there is a perturbative Hartle's 
method which is based on the assumption that rotating massive body is no longer 
spherically symmetric. It is destorted, thus expanding the metric functions through 
second order in the stars rotational velocity Q one can obtain the following form of the 
perturbed metric 

ds 2 = g^dx"dx v = -e 2l/(r ' 9 ' n) c 2 dt 2 + e 2 ^ e ' n) (d(f) - w(r, 9, tycdt) 2 
+ e 2 ^ e ' n W + e 2X ^ e ' n W 2 + 0(fi 3 ), 
where metric functions in this perturbed line element are given by 

e 2u(r, 0, «) = e 2*(r) (1 + 3^ Q) + ^ fi )P 2 ( COS 0))), 
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Figure 8. The mass-radius diagram for the protoneutron star for constant Yp = 
0.1776. 
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Figure 9. The mass-radius diagram for the protoneutron star with T — 30 MeV 
(Y L = 0.4, Y P = 0.1776) and T = MeV( the proton fraction is constrained by (3 
decay). 

e 2^(r, e, n) = r 2 sin 2Q(i + 2 (v 2 (r, Q) - h 2 {r, fi))P a (cos0)), 
e 2 M (r, 6, a) = r 2 (1 + 2 ^ 2 ^ fl) _ ^ fi ))p 2 ( co ^)) ) 

e 2A(r, e, n) = 2A(r) A , 2G m (r, n)+m 2 (r, f])P 2 (cosfl) \ 

where $(r) and A(r) are the metric functions of a spherically symmetric star, P 2 the 
Legendre polynomial of order 2, mo,m 2 , /io, ^2 an d v 2 are all functions of r and fl 
|19| , ||. They are calculated from Einstein's field equations and given as solutions of 
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Figure 10. The time evolution (t = — 15s) of the protoneutron star (with 
MB = IMq) in accordance with paper Toki ||25| . 




Figure 11. The mass M dependence on the central density p c for Yl = 0.4. 



Hartle's stellar structure equations, u has the same meaning as in the nonperturbative 
line element fl35l) . The metric functions are determined with the use of the Einstein 
equations 

G = R^ u — —g^ u R = kT^, (37) 
where is the stress-energy tensor given in the perfect fluid form 

Tpv = (e + P)u^u v + g^P (38) 
and is the unite four-velocity satisfying the following condition 
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Figure 12. The mass M dependence on the central density p c for Yp = 0.1776. 



where e is the total energy density, P the pressure. Each metric function, namely v, p 
and io are functions of the radial coordinate r, polar angle 9 and also the stars angular 
velocity Q. In the case of rotating stars there is an additional dependence of the metric 
on the polar angle and the frame dragging frequency io. The latter leads to the existence 
of the non-diagonal term g ttp in the metric tensor. The assumption of uniform rotation 
means that the value of Q is constant throughout the star. For uniformly rotating 
bodies there is a relation between components of the four- velocity vector = VLu 1 . 
The nonzero components of the four-velocity vector of the matter are of the form 



where 



u 



CO 



n 



10 



co 2 e 2 ^) 



alt 



(39) 



(40) 



The absolute limit on stable neutron star rotation is the Kepler frequency Qk- It 
determines the frequency at which the mass shedding at the stellar equator sets 
in. The result of the work of Haensel and Zdunik 261 shows that the value of the 



Kepler frequency can be estimated knowing the value of the mass and radius of the 
corresponding nonrotating star and an empirical relation was given 

tt K w C Hz ^J(M s /M & )(R s /10km) 3 = (0.63 - 0.67) x fi c (41) 

where, Ch z — 7700s -1 and fl c is the Newtonian value and is equal 

Q c = y/MjjF, (42) 

the index s indicates that these values refer to the spherical configuration. As a 
consequence of the perturbative method for the angular velocity of the local inertial 
frame appears 

d 
dr 



4-/ ,du(r) 
r J{r)— 7— 
dr 



+ 4r 3 



dj{r)_, , 
Lo(r) 



dr 



0. 



(43) 
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Figure 13. The moment of inertia /, (in units of Mq km 2 ) as a function of the central 
density for the different temperatures. The lepton number equals 0.4. 



where cu = Q — uj and j(r) = e — ^^^-)2 and the boundary conditions are 

such that 57(0) = uJ c and = 0. The angular velocity dependence on the angular 

velocity in the center is stated as 

\ / r=R s 

The solution of equation ( (4"3|) allows us to determine the star's momentum of inertia 
m_S* f- drri £W + P(r) gr) e _ t(r . 



n ~" 3c 2 7o (i - 2Gm(r!a n 

The f2 c is angular velocity in the approximation the amount needed to produce shedding 
of mass at the star's equator, so that this method of computation is not actually valid 
for this large a value of fl. 



4. The numerical results. 



In this paper the behaviour of a protoneutron star in the RMF approach was examined. 
In Fig.l the pressure P as a function of the energy density is shown for the TM1 
parameter set in the case of finite temperature. The relevant temperature range is 
between 1 MeV up to 75 MeV. Figures 2 and 3 display the lepton fraction Yl = 
(n e + n u )/ns and the entropy as a function of central energy density. For different 
temperatures the proton fraction Yp is fixed and equal Yp = 0.1776. The pressure 
as the function of the star radius is presented in Fig.4. The total pressure of the 
protoneutron star being the sum of fermion and meson parts is denoted by Pf u ii- 
Particular contributions P u , P p , P n , P e , P ou t, P mark the pressure coming from 
neutrinos, protons, neutrons, electrons, thermal plasma and pressure without neutrinos, 
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Figure 14. The moment of inertia /, (in units of Mq km?) as a function of the central 
density for the different temperatures. The proton number equals 0.1776. 
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Figure 15. The mass-radius relation for the protoneutron star for different value of 
f2. The temperature T and the lepton number Yl are fixed and equal 50 MeV and 
0.4, respectively. 

respectively. Neutrons pressure is the most significant component of the total pressure 
whereas contributions coming from neutrinos and protons are nearly the same in the 
star core. As the radius increases the neutrino pressure starts to prevail over the 
remaining pressures. Figs. 5 and 6 depict the lepton number Yl and the entropy per 
baryon s as the function of the radius. The temperature is fixed and equals 50 MeV, 
p c = 1.8 x 10 15 g/cm 3 . For the considerable value of the star radius the lepton number 
and the entropy are almost constant whereas in the outer layers of the protoneutron star 
both of them grow very steeply. The mass-radius relations for protoneutron stars for 
different temperature cases and constant lepton Yl and proton Yp number is presented 
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Figure 16. The mass-radius relation for the protoneutron star for different value of 
f2. The temperature T and the proton fraction Yp are fixed and equal 50 MeV and 
0.1776, respectively. 




Figure 17. The mass-central density relation for the protoneutron star for different 
value of fl. The temperature T and the lepton number Yl are fixed and equal 50 MeV 
and 0.4, respectively. 



in the next two figures. The stable and unstable areas of the stellar configuration 
(the dotted line) are presented. The possible stellar evolution tracks with the constant 
baryon number are pointed by arrows. Considering configurations obtained for different 
temperatures one can see the possible evolution tracks. They are presented in Table [| 
and Table || At the initial time t = s the protoneutron star is characterized by the 
barion mass of the value of 1 M & and the proton fraction Yp = 0.2, the temperature 
equals 30 MeV. In accordance with the radial profiles of temperature present in paper 
p5fl after 15 s the temperature and the proton fraction drop to 10 MeV and 1/9, 
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Figure 18. The mass-central density relation for the protoneutron star for different 
value of fi. The temperature T and the proton fraction Yp are fixed and equal 50 
MeV and 0.1776, respectively. 



T[MeV] 


R [km] 


M[M ] 


50 


20.8981 


1.9751 


30 


79.9505 
40.6529 


3.2871 
1.9139 


10 


40.5229 


1.6566 


1 


20.2331 


1.6328 



Table 2. The evolution tracks for Y L = 0.4. 



T[MeV] 


R [km] 


M[M e ) 


50 


22.1870 


1.5560 


30 


84.6644 
32.2714 


3.1236 
1.3731 


10 


25.2496 


1.2360 


1 


15.8983 


1.2232 



Table 3. The evolution tracks for Y P = 0.1776. 



respectively. For comparison the mass-radius dependence for the protoneutron star with 
constant lepton and proton number and in the temperature equal 30 MeV are presented 
in Fig.9. The unstable configurations appear and they are marked by dotted lines. The 
same relation is shown for the star in zero temperature limit and with zero neutrino 
chemical potential (without neutrinos), than the proton concentration is constrained by 
the j3 equilibrium. The mass-central density (p c ) functions for different temperatures 
and fixed lepton number Yl (Fig. 11) and proton fraction Yp (Fig. 12) are shown on 
Figs. 11 and 12. From relations presented on Figs. 7 and 8 one can see the density ranges 
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T [ MeV ] 


R km 


M[M Q ] 


Mr [ Mn 1 


I [ M Q km 2 ] 


1 


14.3886 


2.1312 


2.3057 


130.4967 


10 


18.3521 


2.135 


2.3051 


131.3594 


30 


80.0136 
19.6225 


3.3723 
2.1694 


1.7768 
2.2806 


123.7578 
146.7724 


50 


16.4648 


2.1659 


2.1226 


152.4982 


75 


12.5894 


1.9839 


1.7173 


116.9361 



Table 4. The mass, the radius and the moment of inertia for Yl = 0.4. 



T[MeV] 


R [km] 


R [km] 


M[M Q ] 


M[M Q ] 


J[M Q km 2 s' 1 ] 


1 


14.3886 


15.6158 


2.1312 


2.4281 


1271590 


10 


18.3521 


20.9522 


2.135 


2.2798 


889378 


30 


80.0136 
19.6225 


92.0202 
22.5576 


3.3722 
2.1694 


3.7851 
2.3129 


6658570 
906022 


50 


16.4648 


18.2002 


2.1659 


2.4031 


1223790 


75 


12.5894 


13.4091 


1.9839 


2.3060 


1343250 



Table 5. The mass, the radius and the angular momentum for Yl = 0.4. 

for which the stars are stable. Using the results of Figs. 7 and 8 the density ranges of 
unstable protoneutron stars configurations can be estimated. The moment of inertia (I) 
of a protoneutron star versus the central density (p c ) for different values of temperature 
are presented. Figures 15 and 16 show the subsequence mass-radius relations of a hot 
protoneutron star (T = 50 MeV) obtained for different angular velocities. The lepton 
and proton number are constant. The greatest values of the angular velocity the more 
massive and bigger are the stars. Figures 17 and 18 compare the mass-central density 
relations for protoneutron stars with fixed values of lepton Yl and proton Yp number. 
The temperature equals 50 MeV obtained for different angular velocities. In Fig. 18 the 
density range of stable rotating protoneutron star is visable. The main numerical results 
is collected in the Table |] and the Table [5], where Yl = 0.4 and p c = 5 10 14 g/cm 3 . 

5. The conclusion. 

The main aim of this paper was to study the protoneutron star parameters especially 
the masses and radii as the most sensitive ones to the form the equation of state. The 
employed form of EOS in the TM1 set which was extended to the finite temperature 
cases. The considered model comprises not only nucleons and electrons, which are 
necessary for the f3 stable matter, but the neutrinos as well. The presence of neutrinos 
is characteristic for hot young protoneutron stars. Neutrinos give not negligible 
contribution to the pressure and energy density. As the radius of protoneutron stars 
increases the neutrino pressure starts to prevail over the pressure coming from the 
remaining components. Constructing the mass-radius relations for protoneutron stars 
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the stable and unstable areas appeared. This indicate that there are several possible 
stellar evolution tracks for the configuration with constant baryon number. These 
evolutionary tracks are strictly connected with the decreasing star temperature and 
thus with the decreasing neutrino chemical potential and lepton number. The final 
result is the cold deleptonized (the neutrino chemical potential equals zero) object. The 
influence of rotation on the protoneutron star parameters is significant. In this paper 
the case of slowly rotating protoneutron star is considered. The obtained unstable 
configuration for the fixed temperature (T = 30 MeV) and lepton number (Yl = 0.4) 
occurring in the density range (1.4 — 6.5 x 10 14 g/cm 3 ) are relevant to rapid change in the 
moment of inertia of the star (Fig. 13 and Fig. 14). For the rotating star the mass-radius 
relation depends on the star angular velocity thus the values of Q affects the instability 
areas. For the employed EOS there exists the ranges of masses where hot, lepton rich 
protoneutron star can be stabilized against gravity but loosing leptons these objects 
became gravitationally unstable. The instability strictly connected with the 

presence of neutrinos and the thermal pressure effects. 
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